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The equation for the gluon propagator in the approach of Baker-Ball- 
Zachariasen is considered. The possibihty of non-integer power infrared be- 
haviour is studied, D{q) ~ (i?^)"^, — ^ 0. It is shown that the characteristic 
equation for the exponent has no solutions at — 1 < c < 3. The approx- 
imations made to obtain the closed integral equation are analysed and the 
conclusion on the infrared behaviour of the gluon propagator D{q) ~ 

is made when the transverse part of the triple gluon vertex is taken 
into account. 
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The study of the infrared behaviour of Green's functions without any 
doubt is an important problem in QCD. We mean first of all the infrared 
behaviour of the gluon propagator. This behaviour is considered to have a 
decisive role in the problem of confinement, calculation of vacuum expecta- 
tion values of gluon and quark fields, the description of spontaneous chiral 
symmetry breaking, etc. [jl], 0. The renormalization group method allows 
one to obtain the asymptotic behaviour of the gluon propagator at large 
momentum ||^. As we consider the infrared behaviour, then the approxi- 
mate solution of the exact Schwinger - Dyson equations for Green's functions 
seems to be the most adequate method of nonperturbative nonlattice study. 
At present the commonly accepted result on the character of the infrared 
behaviour of the gluon propagator is absent in the literature. 

Nonuniqueness of the solutions can be conditioned by nonlinearity of 
the corresponding integral equations as well as the difference of the trun- 
cation methods of the equations for Green's functions. As it is shown in 
Refs. Q, when finding nonperturbative solutions (and the infrared region 
is completely nonperturbative one) the problem of supplementary bound- 
ary conditions is of great importance. In Refs. [5-9] the gluon propagator 
behaviour D{q) ~ was obtained. This behaviour gives a 

linear confining qq potential at large distances in the Born approximation. 
Moreover this behaviour serves as a basis for physically attractive picture 
where the QCD vacuum is considered as chromomagnetic superconducting 
medium [|lOl. In this paper we consider the equation whose investigation 
allowed the authors of |11| to come to the conclusion on the existence of the 
solution having soft power behaviour at the origin and study the general case 
of power infrared behaviour with non-integer exponents. Then we compare 
the obtained results with those of dealing with the search for the 
field equations solutions for the propagator with non-integer power infrared 
behaviour. Then we analyse the approximations made when obtaining the 
closed equation for the propagator and discuss the results. 

The gluon propagator in Baker - Ball - Zachariasen approach [^] is defined 
by only one scalar function Z{q^)^ 

D,M) = Z{q^)Df}{q\ Df}{q) = -{q')-'^,M, (1) 
If we assume that the longitudinal part of triple gluon function is dominant 
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and neglect the transverse part then for Z{q^) the following equation can be 
obtained [g 0: 



{Z-\q^)-l){l-y)=2-'Z-\q') j 



xSap(A;') [z{k'^){k'n)p{k,q)qM + k)p - 
-Z{k^){k'r^)p{k'q)q,{q + k\ - Z{q^) (Z{k''){kr^) - Zik'){k'v)) g^,~ 



Ziq')pik, k') iikr^) - ik'v)) m')9.p - KK) 



(3) 



In Eqs. (|T|) - (H) 7]^ is gauge vector, k' = q — k, dk = ig'^fi'^~^Nd"k/{2-K)^, n 
is space-time dimension, p{qi,qj) = {Z{qf) — Z{q'j))/{qf — g|). Provided the 
condition [qr]) = is imposed, after a Wick rotation, approximate angular 
integration and the corresponding UV - subtractions, the equation for 
renormalized function Z^i^q^) = Z{—q^)/Z{—p?) has been obtained [0, p^ 
is renormalization point. This equation has the following form: 

1 3aJ/i2) / T2 \ 



ZR{q' 



TT 



ZR{q' 



y /u. 00 00 

T^=j Fr{q\y)dy - J F,{p^ ,y)dy + J F^iq"" ,y)dy - J F2{p\y)dy, 

g2 ^2 

00 00 

T^ = j F,{q\y)dy - J F,{p' ,y)dy + J F,{q\y)dy - J F,{p' ,y)dy.{4) 

q2 ^2 

The functions Fi 2{x,y) in Eq. are linear in Zji{x), Zji{y), Zji{x + y), and 
Fi,^{x,y) are quadratic. As a result of numerical study |]ll| the following 
parameterization was suggested to be used in phenomenological applications 



D{q' 



ZrW_ 
g2 



Ci(g7m2)7i + c2(gVm2)72 + L{q^ / w?)ln{\iq^ / m? + A2) ' 



(5) 

It is essential that 71 ~ 0, 22 and 72 — 0, 86 are non- integer parameters. In 

(6) 



Ref. |T5[ the possibility of the infrared behaviour of the form 

Zr{x) = {x/ p^^ [a^ + aix/ p^ + a2{x/ p^Y + 



was considered for the values < c < 1 corresponding to the soft singular 



behaviour of the propagator. In Ref. [16| the wider interval — 1 < c < 3 



was studied, and in Ref. the generahzation was made for the case of 
superposition of the terms with power infrared behaviour, 

N 
n=l 

where functions k = 0,1, N are supposed to be analytic at the 

origin and non-integer parameters numbered as Ci < C2 < ... < cn 

(ci = c — 1, /i(0) = ctg ^). Equation (||) was shown to have no solutions 
with infrared behaviour of form (P), (^), if — 1 < c < 3 and c is non-integer 
(non-half-integer). When doing this, the technique of extracting power terms 
(~ q"^ 0) of the integral equation under consideration defined only 

by the infrared behaviour of the propagator was developed. In the infrared 
region Eq. (^) takes the form: 



- a,{q'/^?^^/\{c)+o{{q'/^^^^^)}, (8) 

where A(c) is characteristic function, P{q^), Q{q^) are some regular at = 
functions which depend on the behaviour of Zulq"^) at all values of the 
argument. At ci > the term with the behaviour (g^)"'^^ is leading at — 0. 
So, for the exponent we have the characteristic equation A(c) = 0. At Ci < 
the term (g^)"^^ in eq. (|^) defined by the infrared region is not leading. 
Nevertheless, the consideration of this case [|l^, [T^ shows that for the infrared 
asymptotics (|) or (0) to be consistent with Eq. (|) it is necessary that 
A(c) = 0. It turns out that the function A(c) at c = -1,-1/2,0,1/2,1,3 
has simple poles and in the intermediate points of the interval [—1,3] there 
are no roots of the equation A(c) = 0. The explicit form of the characteristic 
function is 

23 1 187 1 7 1 5 1 

A(c) = \ h 

^ ^ 24 c + 1 96 2c + 1 24 c 12 2c - 1 

2 1 3 1 

+ hA(c), 

3c- 1 4c-3 ^ 
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1 

A(c) = —{J dt[i^ + ty{28t-^+^' + 28f-' - 9t-^+' - 16(1 + ^-3+2=) x 



X (1 - t^-')/{i - 1)) - i2r3+2- + 4(1 + 3c)r'+2'= + 

+ 2(16 - 5c - 3c2)t"^+^' + 2(c^ + 4c2 - 13c + 16)t'^^ + 

+ 1 /6(192 - 122c + Slc^ - 22c^ - 3c^)t^+2^ - 16t^~^ + 

+ 4(4c - 15)t^-' - 7t-^+' + (7c - 23)t-^+" - (Tc' - 

- 39c + 64)t72] - 237(2 - c) + 417(3 - c) - 117(4 - c) + 

+ 7(5 - c) + - + 63/4 - 17c/6 + 7cVl2 + cV4}, (9) 

where 7(2;) = (2^ — l)/x. In Ref. [JT^I the approximate angular integra- 
tion in eq. (^ was not carried out and dimensional regularization was used 
for which the ultraviolet subtractions are made automatically. Denote by 



Arfjm(c) the characteristic function of Ref. |jT2[ divided by — 2r(c) to pass to 



the normalization of the present paper. It looks like: 



2r(c + i) 

_ (c^-6c+15)r(c-2)r^(3-c) 

2r(c + i)r(6-2c) ■ ^ ^ 

At c < 2 the behaviour of the function Adim{c) is similar to A(c) but at c > 2 
the behaviour of Adim{c) qualitatively differs IjT^ from that of A(c) and at 
c ~ 2, 537 it vanishes. This value corresponds to nearly oscillatory potential 
V{r) ~ r^'^"^. However, we are not confident that this solution is consistent 
for nonzero values of gauge parameter y = 7"^ = {priY /p^rf . 



In Ref. [|T3| the equation for the gluon propagator in simple loops ap- 
proximation was investigated and self-consistency equation for non-integer 
exponents was obtained. Solutions of this equation were found but it re- 
mains unclear whether its solutions correspond to the infrared asymptotics 
or to the ultraviolet one. 

Let us consider the assumptions made when obtaining the basic equa- 
tion (^ in some detail. The Schwinger - Dyson equation for the gluon 
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propagator in the axial gauge has the following structure: 

PM - pS){p) = ^Ap) = KM + KM, (11) 

where IIj^^, and 11^^ are one-loop and two-loop terms of polarization operator. 
The inverse propagator P^uij)) depending on two scalar functions and free 
inverse propagator Pj^J can be represented by 



where 

= 9tiu - Pf.Pu/p'^, L^=Pf,- r]^p^/ [pr]). 
The tensor structure of one and two - loop terms is 0]: 

KM = Kt;:. + KL,LJp' + UlN,LJipr^), (12) 

KAp) = KkM + Kn,l,/{pv), (13) 



where 



K,. = s^.-7r;, = (i-7)r;. + i^MW/, 

' tJ,u 



= Vf^-Pt^v'^/ipr]), TJl^ = 9]^v -ri],r}u/r{^- (14) 



Note that the presence of nonsymmetric term in Eq. ([T^ ) makes it incon- 
sistent to consider the tensor equation without two-loop term. Let us see 
further what are the consequences of the assumption on dominance of one 
tensor structure in the propagator, D^^{p) = Z{p)D^^J{p). In this case 
Pa = —Z^^{p)p'^, Pb = and tensor Eq. (|ll]) gives Ilf = —11^, so we ob- 
tain the following tensor equation: 

p,M-p!.'J{p) = (K-{i-i)K)T^.= 



Consider the contribution of the tadpole term. We have 



'^K>T'- (15) 



7]^ (27r) 



-/N j<rk^{n-2 + ^). (16) 
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It does not depend on the external momentum and it is some constant, may 
be infinite. Consider the assumption on independence of Z of gauge parame- 
ter, Z{])) = Z{p'^). This assumption has been shown |TB| to be inconsistent in 
the framework of the BBZ approach and this approach should be modi- 
fied. Consider what the consequences of the assumption for the tadpole term 
are. The angular integral equals 



dVLr 



k'^-rf 
{kr]y 



dx 



-(n-2)fi„, (17) 



where f2„ = 27r"/^/r(n/2), and the singularity at the origin is understood in 
any of the senses: 



1 



PV 



a;2' (x + i0)2' {x-i^y 



X 



So, for reasonable axial singularity prescriptions and sensible regularizations 
(e.g., dimensional regularization, symmetric integration) the tadpole term 
vanishes for Z{p) = Z{p'^). For Z{k'^) |fc2_^oo~ In^k"^ the divergent term in 
Eq. (^ of the form ~ A^Zn"A^ should be subtracted by the corresponding 
mass counterterm. The characteristic equations have no solutions at c = 2 
{D ~ l/(g^)^). In this case the infrared contribution is a constant and it 
could be considered as a part of mass renormalization constant |jl8|. Before 
doing this let us look more attentively at the condition {pi]) = imposed to 
make the problem tractable. As it was noted in the Poincare-Bertrand 
singular terms should be taken into account when using the decomposition 
formula for the axial denominators. For generalized c-prescription pO| 



{kr]) 



{kr]) + iO 



c5{{krj)) 



(18) 



the decomposition formula has the form: 



1 



{kri) 



1 



{p - k, 7]) 



(pv) 



{kr]) 



+ 



{p - k, f]) 



+ 



+ C{2TII-C)5{{k7]))5{{p7])). 



(19) 



At c = 0, 27^^ the Poincare-Bertrand terms are absent and in the case c = m 
corresponding to principal value prescription, one should account for them. 
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Consider whether these terms prevent from going to the hmit {prj) = in the 
integral equation (Q). It can be seen that in product of the axial singularities 
we can not go to the limit (prj) = 0, 



PV-^PV- ^ 



{kr]) {p-k,ri) 



^ -PVjl^ - ii^v)) ^ iiPV)) ■ (20) 



{kriY 



It can be shown that for Z{q^) = —M'^/q^ the additional contribution in 
Eq. (^) has the form: 



Apb(p) = g^M^N{pyfiy./^6{{pr]))x 

X {y-^/^+n0i(6) + #2(e) + Oiy')] + xi(e) + + Oiy')}, (21) 

and 01(e) = 0, Xi(e) |e-^o~ e, e = (n - 4)/2. At e = 

ApB{p) = 6{{pv))J2ck{pv)' = 0, 

k=l 

and the Poincare-Bertrand terms do not prevent from going to the limit 
{prj) — > 0. It is reasonable to assume that this is true for the less singular 
behaviour of Z{k^) at zero. 

Consider the results for the case [prj) ^ 0. For Zlq"^) = —M'^/q^ it was 
obtained in Refs. [H, 0, pll] that, 



nj.^ = -^^[3^ - 3 - yin{^y) - ny)/% (22) 



where 

F{y) =y r lnity)il - t)-'/^il - ty)-'dt. 



If one keeps the dependence on e, then aXy ^ the result can be represented 
as a series of the form 



nj.^ = -^^[7 - W + a,{e)y + h{e)y'^^ + ■■■]. (23) 



At e < it is impossible to go to the limit y 0, and at e > one can do 
it with the result 7 in brackets of Eq. (p3|) . If one puts e = at first, then 
—3 in brackets of Eq. (|23| ) at y = is obtained. Thus the limits y ^ and 



8 



e — s> do not commute and the point {prj) = is singular for Eq. (j^) at least 
for ~ l/q"^- Having the dependence on the gauge parameter y in (|2^) , 

it is impossible to solve Eq. asymptotically by mass subtraction. 

Let us verify the correctness of the assumption on dominance of longitu- 
dinal part of the vertex function. The direct calculation ^ shows that the 
transverse part of triple gluon vertex introduced in ref. [0 of the form 

^'^Ixi.uiP^ r) = ^M-'^{9xi,[Pu{qr) - qu{pr)] + 

+ 9t^u[qxipr) - rxipq)] + gvx[r,,{pq) - pM^)] + rxp^^qy - qxr^^Pu} (24) 

gives nonzero contribution which completely cancels y-dependent right hand 
side of Eq. (|^) at ,^ = 1. So, taking into account the transverse part of 
triple gluon vertex (which is not fixed by Slavnov - Taylor identities) it turns 
out possible to satisfy asymptotically the Schwinger - Dyson equation for 
the propagator with infrared behaviour D{q) ~ l/{q^Y,q^ 0. Note that 
the transverse part of vertex function ( p4D naturally meets the scaling prop- 
erty pri . It seems unlikely that there are other ways to cancel the leading at 



^ term of the Schwinger - Dyson equation which has the complicated 
dependence on the gauge parameter. 

We conclude that the power type infrared behaviour (|^), (|^) of the gluon 
propagator is inconsistent, whereas the behaviour D{q) ~ I / [q"^)"^ , q"^ — >■ is 
consistent. Despite of the differences in the argumentation, the same conclu- 
sion was made by the authors of the recent paper 
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